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The spatial organization and time dependence of Jupiter’s temperatures near 250-millibar 
pressure were measured through a jovian year by imaging thermal emission at 18 mi- 
crometers. The temperature field is influenced by seasonal radiative forcing, and its banded 
organization is closely correlated with the visible cloud field. Evidence was found for a 
quasi-periodic oscillation of temperatures in the Equatorial Zone, a correlation between 
tropospheric and stratospheric waves in the North Equatorial Belt, and slowly moving 
thermal features in the North and South Equatorial Belts. There appears to be no common 
relation between temporal changes of temperature and changes in the visual albedo of the 
various axisymmetric bands. 


Planetary atmospheres are characterized in 
a fundamental way by their temperature 
structures that govern dynamical forces, 
reflect the influences of radiative and other 
sources of heating and cooling, and estab- 
lish the context in which gaseous chemistry 
and cloud microphysics take place. We 
describe the first results of a study of the 
spatial morphology of Jupiter’s temperatures 
through a jovian year. This work probes a 
region where energy transfer is influenced 
by convection as well as radiation, deeper 
in the atmosphere than addressed in our 
study of stratospheric temperatures ( 1 ) . The 
global and long-term nature of our results is 
a major extension of previous studies of 
Jupiter’s tropospheric temperatures and dy- 
namics, which have largely been derived 
from short-duration spacecraft studies (2). 
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We observed Jupiter between 1980 and 
1993 at the NASA/Infrared Telescope Fa- 
cility on Mauna at Kea, Hawaii (Table 1). 
Tropospheric temperatures were deter- 
mined from thermal infrared emission ob- 
served by means of a facility photometer 
with a discrete filter centered at 18.20 pan 
(549.4 cm -1 ) (3). Near 18 |im, Jupiter’s 
atmospheric opacity is dominated by the 
collision-induced absorption of H 2 (4) and 
is minimally influenced by clouds. As H 2 is 
uniformly mixed throughout the atmo- 
sphere, the substantial temporal and spatial 
variations of the thermal emission observed 
with this filter (Fig. 1) indicate variations of 
atmospheric temperature. 

Central meridian scans obtained in 1980 
through 1982 were sampled at 1 arc sec 
intervals. Partial coverage of Jupiter’s disk 
followed in raster scans made in 1983 and 
full coverage thereafter, sampling at 1 arc 
sec intervals in right ascension and declina- 
tion. The diameter of the circular aperture 
used in 1980 and 1981 was 3.92 arc sec. For 
all later observations, the diameter was 
1.96 arc sec, close to the diffraction limit, 
corresponding to 5° or 6600 km at the 
equator. The assignment of latitude, longi- 
tude, and emission angle to each pixel of an 
image required accounting for the drift of 
the telescope from Jupiter’s actual motion 
in the sky and its rotation during the 
approximately 20 to 45 minutes required to 
scan the disk. Spatial uncertainties are on 
the order of 1 arc sec or less. Our study 
includes only locations on the disk with 
emission angles of 56° or less, as positions 


closer to the planetary limb suffered percep- 
tibly from blending with cold space. 

Just as with our 7.8-|Jim images (1), we 
derived absolute radiances by normalizing the 
mean value observed in the inner 50 percent 
of a strip along Jupiter’s central meridian to 
the average radiance observed by the Voyager 
1 IRIS central meridian maps in a simulated 
response to the 18- pan filter. The net effect 
minimized temperature changes with time for 
the largest number of points on the planet, 
rendering the data insensitive to changes of 
the global mean temperature, but preserving 
relative differences. The best images have a 
signal-to-noise ratio of about 400 (equivalent 
to an uncertainty of just under 0.05 K in 
brightness temperature). Noisier images re- 
sulted from poor seeing or from local humidity 
variations which map the variable signal onto 
the observed disk as a false spatial variation. 
The signal-to-noise ratio for most of the im- 
ages is closer to 100 (0.15 K in brightness 
temperature). We did not use images with 
signal-to-noise ratios less than 70. 

In order to interpret the observed radi- 
ance in terms of physical temperature, we 
assumed an average vertical temperature 
structure. A smoothed average of Voyager 1 
radio occultation experiment results (5) was 
used in an initial guess of the temperature 
profile; the overlying stratospheric temper- 
atures were cooled, however, (Fig. 2) in 
order to provide the best fit to the radiances 
observed across the disk at all emission 
angles. The resulting profile is close to that 
of the Voyager radio occultation egress 
profile (5). This profile was perturbed by 
the same temperature at all levels among 
the family of profiles illustrated in Fig. 2 
until the observed radiance was matched; 
this provided a one-to-one correspondence 
between a temperature profile and the in- 
tensity at each pixel of an image. While 
temperature perturbations that decrease 
slightly with depth might correspond more 
closely to Voyager results (5), the differ- 
ences between the temperatures derived 
from such perturbations and the simple 
scheme adopted here are not significant. 
For locations observed at more than one 
emission angle, the temperature was de- 
rived from the observation at the smallest 
emission angle, or, in some cases, from the 
least noisy map. Temperatures are reported 
at 250 mbar, near the peak of the contri- 
bution functions for nadir and limb obser- 
vations (Fig. 2). While this simple proce- 
dure is insensitive to possible variations of 
the vertical temperature profile from those 
shown in Fig. 2, such variations are not 
expected, as the different radio occultation 
profiles show strikingly similar morpholo- 
gies near 250 mbar. 

Zonal mean thermal structure. At any 
given latitude, Jupiter exhibits significant 
temporal and spatial variation of its zonal 
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mean (longitudinally averaged) thermal 
structure. Selected profiles of the zonal 
mean 250-mbar temperature between 1978 
and 1993, together with profiles of the 
vertical shear of the zonal mean wind are 
shown in Fig. 3. Most of these profiles were 
constructed from separate images (Table 1) 
acquired within 3 days and covering a large 
longitude range, although some epochs 
have large gaps in temporal and longitudi- 
nal coverage. The temperature profile for 
1979 was derived by simulating the 18-|xm 
filter response to Voyager IRIS north-south 
map sequence spectra. The profile shown 
for 1978 was derived from the 18-|xm cen- 
tral meridian scan of Sinton et al. (6) . The 
shear was derived from the temperature 
profile using the thermal wind equation, 

du R dT 

^ dz aH Scj) 

where u is the eastward zonal mean wind, T 
is the temperature, z is minus scale height 
times log of the pressure and is used as a 
vertical coordinate, 4> is latitude and /, a, 
H, and R are the Coriolis parameter, plan- 
etary radius, pressure scale height, and gas 
constant, respectively. The geostrophic ap- 
proximation on which Eq. 1 is based breaks 
down near the equator where f becomes 
small, but because Jupiter is a rapidly rotat- 


Table 1. Observational summary. 


Dates (UT) 

Scans (no.) 

1978 February, from (6) 

2 

1980 January 

2 

1981 March 

2 

1982 February 

2 

Dates (UT) 

Maps (no.) 

1983 July 15, 17 

2 

1984 June 6, 8 

2 

1 984 July 23 

2 

1984 August 9, 10 

2 

1984 September 10 

1 

1985 May 21 

1 

1985 September 5 

1 

1986 June 6 

1 

1986 July 9 

1 

1987 June 7, 8 

2 

1987 October 25, 26, 27 

4 

1987 December 30 

2 

1988 September 25, 26, 28 

4 

1988 November 30 


December 1 , 2 

3 

1989 December 15, 16 

5 

1 990 January 9 

1 

1990 February 11 

1 

1990 April 6 

1 

1991 February 9, 10 

3 

1991 April 2 

1 

1991 May 5 

1 

1992 October 26, 27 

2 

1993 March 1 

2 

1993 March 30, 31 

1 

1993 April 27 

3 


ing planet with zonal winds that are steady 
over one rotation period, we expect the 
equation to yield accurate wind shears to 
within a few degrees of the equator. We 
estimated the zonal wind shear at 0° lati- 
tude from the formula Hdu/ dz = — [R/ 
(2fta)]d 2 T/d(|) 2 (the planetary rotation pe- 
riod is 2t t/H) , which strictly is valid only in 
the case when the temperature profile is 
symmetric about the equator. 

The time-averaged temperature and ver- 
tical wind shear profiles reinforce the inter- 
pretation of the IRIS data (7). The North 
and South Equatorial Belts (NEB and SEB) 
near ±12° are some 2 to 3 K warmer than 
the Equatorial Zone (EZ) and mid-latitude 
zones. As found by Gierasch et al. (7), the 
vertical shear of zonal wind tends to be 
anti-correlated with the zonal jets at the 


cloud-top level, implying a decay of zonal 
wind with height above the clouds. 

Most of the temporal behavior evident 
in Fig. 3 consists of (i) a seasonal variation, 
(ii) a quasi-periodic oscillation at the equa- 
tor, and (iii) strong variation of tempera- 
tures and vertical wind shears in the lati- 
tude band 20 to 27°N. The time variation 
of temperatures at eight specific latitudes 
from 1978 to 1993 is shown in Fig. 4- 
Jupiter’s orbital period is approximately 
11.8 years; the last northern summer sol- 
stice occurred in 1988.7. The influence of 
seasonally varying insolation is most evi- 
dent at high latitudes, where more or less 
steady warming is evident between 1984 
and 1989 during northern hemisphere 
spring and summer, while southern lati- 
tudes cooled during the same period. Sea- 



Fig. 1. Raster-scan im- 
ages of Jupiter, repre- 
senting the best of the 
images available for 
each apparition since 
1984, identified by the 
epoch of observation. 
Detailed observation 
dates are given in Ta- 
ble 1 . The warmest re- 
gions in each image 
are shown in white and 
the coolest in dark red. 


Fig. 2. The end members 
and center of a family of tem- 
perature profiles used to 
generate a table relating out- 
going intensity as a function 
of emission angle with the 
temperature at a given atmo- 
spheric level. These have a 
cooler stratosphere than a 
nominal profile, shown by 
the dashed curve, which 
represents a mean of Voyag- 
er 1 radio occultation exper- 
iment results (5). Also shown 
as dotted curves are the lev- 
el by level contributions to 
the outgoing thermal intensi- 
ty as a function of height for 
the baseline temperature 
profile for the 18 -|i.m filter at 

the nadir (lower curve) and at an emission angle of 56° (upper curve). The contribution functions 
have arbitrarily normalized amplitudes but are correct relative to one another. 



Temperature (K) 
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sonal maxima and minima appear within 2 
years after the solstice, which is roughly con- 
sistent with the radiative time constant of ~6 
X 10 7 s (8) which implies a lag between the 
insolation and the thermal response of —1.5 
year. The amplitude of this seasonal signal is 
roughly consistent with seasonal models (9), 
but it is surprising that the seasonal ampli- 
tudes in the southern hemisphere are higher 
than their northern counterparts. Summer 
solstice in the northern hemisphere occurs 
very near perihelion, suggesting that the 
northern hemisphere should be warmer. The 
annual component of the response to seasonal 
forcing is an important signal in the temporal 
behavior of zonal mean temperatures. Never- 
theless, it is mixed with several of other 
variations, which make it difficult even to see 
the repeating nature of the seasonal cycle for 
some regions. 

Among the variations not related to 
seasonal insolation variations is that of the 
EZ which appears to cool by nearly 3 K 
between 1978 and 1983. Despite the lower 
spatial resolution of the ground-based data 
before 1982, we believe this to correspond 


to a real temperature decrease, since the 
Voyager IRIS experiment measured an 
equatorial temperature decrease of 2 K at 
150 mbar occurring during the 4-month 
interval separating the Voyager 1 and 2 
encounters (7). After 1985, the equatorial 
temperature undergoes approximately two 
cycles of a quasi-periodic oscillation with a 
variable period ranging from 2 to 5 years 
and an amplitude of —0.4 K. We believe 
that this is related to the quasi-periodic 
oscillation discovered in Jupiter’s strato- 
sphere (1), a possible analog of the Earth’s 
quasibiennial oscillation ( 10 ). A compari- 
son of the time variation of the equatorial 
temperature at the 250-mbar and 20-mbar 
levels is shown in Fig. 5 A (upper panel). 
The 20-mbar temperatures correspond to a 
coarse average over a depth of approximate- 
ly two scale heights between 50 and 1 mbar 
( 1 ). The oscillations at the two levels re- 
main roughly 180° out of phase even when 
the period of oscillation grew noticeably 
shorter after 1990. Stratospheric and tropo- 
spheric temperature oscillations were de- 
rived from a two-dimensional (latitude and 


Fig. 3. Zonal mean tempera- 
ture and vertical wind shear 
as a function of planetocen- 
tric latitude and time. Epochs 
are identified by year and 
month {see Table 1). The dot- 
ted curves surrounding each 
(solid) zonal mean tempera- 
ture profile illustrate ±1 <r, 
which includes contributions 
from both the noise and the 
actual longitudinal variations 
of temperature, with the latter 
dominating. The vertical dot- 
ted lines show the location of 
zonal wind shear which mark 
the historical northern and 
southern boundaries of the 
NEB and SEB. The large min- 
imum value of the shear, of- 
ten off the scale, required to 
illustrate more subtle varia- 
tions, is (1982) -42.8, (1983) 
-64.4, (1984) -75.6, (1985) 
-51.5, (1986) -55.8 m 

sec -1 (scale height) -1 . The 
values of the zonal wind 
shears within a few degrees 
of the equator are uncertain 
owing to the breakdown of 
the thermal wind equation 
and the uncertainty in the as- 
signment of latitude. Never- 
theless, at least qualitatively, 
the strong negative shear just 
south of the equator and its 
time variability appear to be 
real. 



-60 -30 0 30 -60 -30 0 30 60 

Planetocentric latitude 


height) dynamical model in which an equa- 
torial temperature oscillation was forced by 
stresses exerted by vertically propagating, 
equatorially trapped waves (Fig. 5B) (11). 
These calculations were very similar to 
those performed in similar models of the 
Earth’s quasibiennial oscillation (12). The 



1980 1985 1990 


Fig. 4. Zonally averaged 250-mbar tempera- 
tures as a function of time for eight different 
latitudes. Temperatures were averaged in the 
2° latitude bins corresponding to 24°N (plane- 
tocentric) for the center of the NTB, 20°N for the 
NTrZ, 12°N for the NEB, 2°N for the EZ, 8°S for 
the northern component of the South Equatorial 
Belt (SEB), and 22°S for the STrZ. In addition, 
temperatures at 46° N and 46° S are shown. The 
position of the NTB, determined from a local 
maximum temperature poleward of the warm 
NEB, appeared to vary with time (see Fig. 3). 
Southern solstice at 1982.76 and northern sol- 
stice at 1988.69 are marked with dotted lines. 
Dates of the observations are given in Table 1 . 
Observations with open symbols denote aver- 
ages over a single map, that is, over only a 
—70° range of longitude. Vertical bars denote 
the standard deviation of the temperature 
mean, arising from both the observational noise 
and the zonal variability of temperature. For 
example, the standard deviations in the EZ are 
much lower than those for the NEB or SEB. 


SCIENCE • VOL. 265 


29 JULY 1994 


627 





ratio of the tropospheric to stratospheric 
oscillation amplitude in the model (0.7) is 
approximately consistent with the data, but 
the relative phasing of the signals at the two 
pressures is different. Comparing the phas- 
ing of the oscillations in the model to the 
data in any detail may not be very mean- 
ingful, however, since the model assumes 
oscillations forced at a single frequency, 
whereas the data indicate a quasi-periodic 
oscillation. More observations and model- 
ing are required to determine whether the 
tropospheric signal is truly related to the 
stratospheric oscillation, but we propose 
such a relation as the best explanation we 
have at present for the behavior of the 
equatorial tropospheric temperature be- 
tween 1985 and 1993. 

The equatorial coolings observed near 
1980 and 1992 (Fig. 4) coincided with 
periods of visible whitening of the EZ (13- 
15). These are consistent with a picture of 
local upwelling during cooling, which 
brings relatively fresh and bright condensed 
NH 3 to a detectable altitude or which 
increases the local amount of saturated NH 3 
gas available for condensation. However, 
the EZ cooling in late 1988 does not corre- 
spond to such a visible change. 

Among the most interesting temporal 
variability seen in the ground-based data 
occurred within a region bounded on the 
south by the northern boundary of the NEB 
(~15°N) and on the north by the northern 
boundary of the NTB (~27°N). In the 
center of this region lies the strong prograde 
jet at 20°N; to north and south lie the 
retrograde jets associated with the northern 
boundaries of the belts. The shape of the 
temperature profiles in this latitude band 
changed significantly with time, especially 
during 1984-1990 (Fig. 3). Accompanying 
the temperature changes in the latitude 
band between the NEB and NTB is a 
substantial reduction of the 250-mbar pos- 
itive vertical zonal wind shear over the 
retrograde jet at the northern boundary of 
the NEB and of the negative vertical shear 
over the strong prograde jet. These vertical 
shears are likely to be underestimated owing 
to the —5° spatial resolution of the mea- 
surements, particularly for the shear over 
the narrow prograde jet (7). If the zonal 
wind at the cloud-top level (—600 to 700 
mbar) is relatively steady (16), the acceler- 
ation at the 250-mbar level of the mean 
zonal wind over the retrograde jet implied 
by the change in the vertical shear amounts 
to at least —3 m sec -1 per terrestrial year. 

The qualitative sense of the acceleration 
over the prograde and retrograde jets suggests 
that they underwent a transition from a state 
of general decay with altitude above the 
cloud-top level to a state in which they were 
nearly constant with height. Presumably, it is 
the result of convergence of the meridional or 

628 


vertical flux of momentum transported by 
turbulent eddies in the region, and hence its 
magnitude and sign may provide clues to the 
nature of the forcing such eddies exert on the 
zonal mean flow. If eddies transporting mo- 
mentum horizontally were absorbed with- 
in the width of the retrogr ade je t (—5°), 
then a horizontal wave stress pu V of order, 
of 0.04 N m -2 is implied (where p is the 
density, u' and v' denote departures of the 
local eastward and northward wind from 
their zonal averages, and the overbar de- 
notes zonal averaging) . This stress is smaller 
by about two orders of magnitude than 
estimates for the rate that kinetic energy is 
fed into the zonal jets at the cloud-top level 
(at about 500 mbar), which are on the 
order of 2 N m~ 2 (17). If instead a vertical 
eddy momentum flux were absorbed over a 
depth of about one scale height, centered at 
the 250- mbar le vel, the implied eddy stress 
would be p u'w' — 10 -4 N m -2 , where w' is 
the departure of the vertical velocity from 
its zonal average. This number is of the 
same order of magnitude as estimates de- 
rived for the vertical stress of equatorially 
trapped waves at Jupiter’s equator (10) and 
of waves at high latitudes (18), and there- 
fore may indicate the general order of mag- 
nitude of wave momentum carried upward 
by vertically propagating waves generated 

Fig. 5. Comparison of the 
time behavior of the equato- 
rial zonal mean temperature 
at the 250-mbar and 20- 
mbar levels. Only the devia- 
tion of temperatures from 
their time average is shown. 

(A) The 250-mbar tempera- 
ture from this study (solid 
line) is compared with the 
20-mbar temperature de- 
rived by Orton et at. ( 1 ) 

(dashed line). (B) Time be- 
havior of zonal mean tem- 
peratures at the 250-mbar 
(solid line) and 20-mbar 
(dashed line) levels, calcu- 
lated with a two-dimensional 
(latitude-height) dynamical 
model (11). Equatorial oscil- 
lations of zonal mean wind 
and temperature in the mod- 
el were forced by a vertically 
propagating, eastward pro- 
pagating kelvin wave and an 
idealized, westward pro- 
pagating "anti-kelvin” wave 
(10). Both modes corre- 
sponded to wave number 1 
and had phase speeds of 
±120 m sec~T The model 
calculation shown produced 
a single oscillation period of 
6 years, while the data ap- 
pear to indicate quasi-peri- 
odic behavior with dominant 
periods ranging between roughly 2 and 6 years. 
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by Jupiter’s convection. The true nature of 
the mechanism that caused the observed 
changes in the vertical wind shears cannot 
be ascertained without detailed observa- 
tions of the eddy winds and temperatures. 

It is also interesting to note changes of 
the visible appearance of these regions 
which were concomitant with the thermal 
changes. The visually dark NTB brightens 
periodically (19, 20) to the same level as the 
nearby North Tropical Zone (NTrZ). It 
brightened in the middle of 1987 and then 
darkened in 1 990 during a major outbreak of 
white and dark spots in its southern compo- 
nent (14, 21-23). The 1987-1990 warming 
(Fig. 4) may be associated with this bright- 
ening (13, 20, 22, 23). The visually dark 
NEB also broadened to the north in early 
1988, and then receded during 1989, leaving 
an array of dark barges and bright ovals (14, 
24, 25). This broadening coincided with a 
pronounced warming in the NEB and the 
latitude labeled as the NTrZ in Fig. 4- The 
South Tropical Zone (STrZ) warming in 
1993 (Figs. 3 and 4) corresponded to the 
most extensive darkening of this zone (15). 
The NEB broadening (or NTrZ darkening) 
and the STrZ darkening are consistent with 
rapid adiabatic warming with subsidence and 
resublimation of condensed (and visibly 
bright) NH 3 ice particles to reveal darker 
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particles below or darker nucleation sites. 
The coincidence of the visible brightening 
of the NTB and its warming or merger with 
the NEB is not so easily understood. How- 
ever, some axisymmetric regions whose vis- 
ible clouds were quite variable exhibited 
little correlated temperature variability. For 
example, there were no thermal variations of 
the STrZ before 1991 during the develop- 
ment of several discrete columns and streaks 
(13, 26 , 27 ). Another surprise was that the 
SEB, which has undergone several dramatic 
albedo changes at all longitudes between 
1989 and 1993 ( 14 , 15 , 25), showed no 
contemporaneous change of its temperature 


structure which is easily separable from the 
seasonal signature. 

Longitudinal structure. Maps of the 250- 
mbar temperature field which are the most 
complete in longitudinal coverage are shown 
in Fig. 6. The most prominent feature in 
these maps is the Great Red Spot (GRS), 
which appears as a very cold oval near 18°S 
in each of the maps, subtending ~ 10° in 
latitude and 18° in longitude, commensurate 
with its visible appearance. The 250-mbar 
temperature of the GRS is 115.3 ± 0.4 K, 
making it the coldest region on the planet 
between 60°S and 60°N latitude. It is about 
1.3 ± 0.5 K colder than the relatively cold 



Fig. 6. Cylindrical projections of temperature profiles derived from selected groups of raster-scan 
images of Jupiter taken over several sequential days between 1984 and 1993 (see Table 1). The 
grid scale is 2° in latitude and longitude. Unusually large excursions from the average near the limbs 
and “seams” between overlapping maps are the result of the inconsistency between the actual 
temperature profile and the assumed form of the temperature profiles (Fig. 2). The temperature 
scale for all the maps is 112.0-121.5 K, except for 1984 August which is 113.0-122.5 K. 


Fig. 7. Plots of stratospheric 
temperature (A, 20-mbar 
pressure) and tropospheric 
temperature (B, 250-mbar 
pressure), for 1988 Septem- 
ber 26-28 and November 
30-December 2, as it varies 
with System III longitude. 
The tropospheric tempera- 
ture profiles correspond to 
13°N latitude (planetocen- 
tric), while the stratospheric 
temperatures correspond to 
22°N. The 1988 November 
temperatures for both strat- 
osphere and troposphere 
are plotted with a longitudi- 
nal offset of -24° in order to 
illustrate the strong correla- 
tion of individual features in 
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axisymmetric STrZ immediately to its south, 
near 22°S, roughly consistent with Voyager 
1 IRIS experiment results at roughly the 
same spatial resolution. 

Longitudinal structure is evident to some 
extent at all latitudes in the maps. It is most 
prominent in the bright NEB, but it is also 
nearly always evident in the SEB and the 
NTB. There are also occasional prominent 
warm features around 35° to 40°S; one of 
these, seen in the 1984 August map, is the 
warmest feature in all of our retrievals. In 
general, to the limits of our resolution, 
thermal eddies at 250 mbar appear to be 
most prevalent at latitudes overlying cloud- 
top retrograde jets. The same result was 
found in analysis of Voyager IRIS data (7). 
The 1989 December map is the only exam- 
ple of identifiable wave structure in the EZ. 
Also, longitudinal wave structure in the 
NEB and SEB is particularly strong between 
1988 and early 1990, and again in 1993. 

Waves observed at 13°N in the NEB 
from September to November 1988 are 
particularly interesting because planetary- 
scale waves were also observed in the strato- 
sphere at 22°N at this time (1). This sug- 
gests the possibility that the stratospheric 
waves might have been the result of upward 
vertical propagation of the tropospheric 
disturbance. The tropospheric 250-mbar 
temperature at 13°N latitude as a function 
of longitude for the dates 1988 September 
26 to 28 and November 30 to December 2 
is shown in Fig. 7. Also shown is the 
stratospheric 20-mbar temperature at 22°N 
for these dates. The November 250-mbar 
temperature is plotted with a longitudinal 
offset of —24° to illustrate the correlation 
between the largest amplitude features for 
the two different times. This longitudinal 
offset corresponds to a zonal group velocity 
for the disturbance of —5.5 m sec -1 mea- 
sured relative to System III, the longitude 
system associated with Jupiter’s radio emis- 
sion and, presumably, the internal rotation. 
The November stratospheric temperature is 
plotted with the same longitudinal offset, 
and we find that relatively large amplitude 
features at —250°, 270°, and 350° longitude 
in the September data appear to have cor- 
responding features in the November data. 
However, the November stratospheric tem- 
peratures are quite noisy, and this precludes 
drawing any firm conclusions regarding the 
correlation of September and November 
features. The relation between features seen 
in the troposphere and stratosphere is also 
shown in Fig. 7. We believe that we can 
identify the relatively warm feature occur- 
ring near 10° longitude in the September 
troposphere with the hot feature near 350° 
longitude in the September stratosphere. 
Similarly, the warm feature at 280° in the 
troposphere probably corresponds to a fea- 
ture at 240° in the stratosphere. The latter 
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identification is more evident in the No- 
vember data. If we accept the premise that 
the tropospheric and stratospheric waves 
correspond to a single disturbance, then the 
identifications made above imply a varia- 
tion of wave phase with altitude, suggesting 
that the disturbance is vertically propagat- 
ing, rather than exponential, in nature. 

Since the thermal waves observed in 
September and November 1988 are com- 
posed predominantly of zonal wavelengths 
of order 40,000 to 70,000 km, and since 
their amplitudes suggest that the Rossby 
number associated with their motion is 
small, it is likely that their propagation is 
governed by the dispersion relation for lin- 
ear, planetary scale Rossby waves. Such 
waves can propagate vertically only when 
their phase speed, measured relative to the 
local zonal mean wind, is opposite in sign to 
the local latitudinal gradient of potential 
vorticity. If the effect of vertical stratifica- 
tion of the atmosphere on the potential 
vorticity q is negligible, then the latitudinal 
gradient of q is given approximately by 

(1 /a)dq/d<$> — (1 /a)df/d<$> — (l/a 2 )d 2, Tr/d<j> 2 

In the retrograde jet at 15°N, which has an 
amplitude of —30 m sec -1 and a width of 
—4° of latitude, (1 /a)dq/d<$> is negative 
(24), and as a consequence, vertical prop- 
agation occurs for waves with phase speeds 
that are eastward relative to the mean zonal 
flow. Our deduction that the thermal waves 
observed at 13°N are vertically propagating 
is consistent with this result; their phase 
speed of —5.5 m sec -1 relative to system III 
is eastward with respect to the zonal wind 
over most of the retrograde jet. (Actually, 
we have deduced the zonal group velocity 
rather than the phase speed from the lon- 
gitudinal offset shown in Fig. 7, but the two 
velocities are very nearly equal for the case 
of Rossby waves with zonal wavelengths 
much longer than the deformation radius, 
which is of order 1000 km or less). Such 
waves are expected to be latitudinally 
trapped in a waveguide within the jet be- 
tween the critical surfaces at which the 
zonal wind matches the phase speed. Nor- 
mally, the wave propagation region would 
also be bounded above by a critical layer 
because the vertical shear is positive (east- 
ward) above the retrograde jet, but during 
times of weak vertical shear, the waves 
might be able to penetrate to high altitudes. 
Such a period of weak vertical shear existed 
during September and November of 1988. 
If this interpretation is correct, we should 
expect to see a correlation between periods 
of high stratospheric wave activity and pe- 
riods of reduced vertical shears over the 
retrograde jets associated with the NEB and 
SEB, a view which is generally consistent 
with the available data but in need of more 
rigorous testing. 


We also examined the data for the pres- 
ence of “slowly moving thermal features”, 
such as those reported in the Voyager IRIS 
data (29). Correlation coefficients were 
computed by normalizing to the longitude 
extent covered by both maps under consid- 
eration, with offsets reported with respect to 
System III. The presence of significant cor- 
relations (in Table 2) is evidence for such 
features. However, with few exceptions, 
they were found only after 1986 and only for 
those latitudes with the highest amplitude 
waves, that is, the NEB and SEB. For the 
SEB, this correlation is always dominated by 
the cold GRS when it is detectable and the 
phase velocity corresponding to these offsets 
divided by the difference in observation time 
is in the range of —4 to — 5 m sec -1 . For the 
NEB, this value of the phase velocity ranges 
between —2 and — 12 m sec -1 . We suspect 
that most of this variability is real. For 1989 
December through early 1990, some wave 
structure is easily detectable at the equator. 
A comparison between the 1990 January and 
1990 February single-image maps shows a 
significant equatorial correlation midway be- 
tween the NEB and the SEB values, although 
we found no significant correlation at the 
equator between the 1989 December maps 
and either 1990 January or Febmary maps. 

Deming et al. (30) detected wave num- 
ber 10 features covering the entire low- 
latitude region of the planet. Our data 
revealed no features of this type, particular- 
ly any covering a broad range of latitudes 
(except possibly in December 1989, as 
shown in Fig. 6), although we saw some 
evidence for slowly moving thermal fea- 
tures, as described above. One possible 
reason is that Deming et al. observed in a 


spectral region which was also sensitive to 
variations of the NH 3 gas and ice distribu- 
tion, as well as to variations of temperature. 
Furthermore, the amplitude of their waves 
was somewhat smaller than the character- 
istic noise level of most of our data. 

Remaining questions. Foremost among 
our remaining questions concerns the rela- 
tion between changes in zonal mean tem- 
perature and visual albedo in major axisym- 
metric regions. The coincidence of a ther- 
mal warming with a decrease of albedo in 
the NTrZ during 1988 is consistent with a 
mechanism where warming of the region 
led to resublimation of visibly bright NH 3 
ice particles, releaving darker particles at 
greater depths or leaving behind dark nu- 
cleation cores. On the other hand, warm- 
ing of the NTB during 1987 to 1990 was 
accompanied by an increase in the visual 
albedo, while albedo changes in the SEB 
were not accompanied by significant tem- 
perature changes, at least at the same time. 
Part of our further inquiry into this matter 
will be to analyze the physical properties of 
clouds from our own infrared database, 
together with a description of the distribu- 
tion of NH 3 gas and tropospheric tempera- 
tures at pressures greater than 250 mbar. 
Another approach to this question is to 
examine the correlation between longitudi- 
nal variability of clouds and temperature, 
which we addressed here only for the axi- 
symmetric regions and the GRS. 

We also have not fully differentiated 
between temporal thermal variations that 
arise from seasonal radiative forcing rather 
than other effects. Modeling the influence of 
seasonal radiative forcing on temperatures in 
some detail will allow us to separate these 


Table 2. Drift rates of 250-mbar temperature waves. 


Epoch 

Planetocentric 

latitude 


Delta 


Drift rate 
(m/s) 

Start 

End 

Longitude 

(°W) 

Time 

(days) 

1987 October 

1987 December 

-10 

-20 ± 

10 

64 

-4.3 ± 2.1 



+ 14 

-12 ± 

6 


-2.6 ± 1.3 

1988 September 

1988 December 

-20 

-20 ± 

4 

66 

-4.4 ± 0.8 



-12 

-16 ± 

2 


-3.4 ± 0.4 



+ 12 

-20 ± 

2 


-4.2 ± 0.4 



+ 18 

-22 ± 

2 


-4.5 ± 0.5 

1989 December 

1990 January 

-10 

-30 ± 

10 

21 

-19.9 ± 6.6 



+ 14 

-4 ± 

4 


-2.6 ± 2.6 



+22 

-10 ± 

4 


-6.2 ± 2.5 

1990 January 

1990 February 

-10 

-30 ± 

2 

33 

-12.7 ± 0.8 



+2 

-28 ± 

2 


-11.8 ± 1.7 



+ 14 

-20 ± 

2 


-8.3 ± 0.8 

1991 February 

1991 April 

+ 14 

-14 ± 

2 

51 

-3.8 ± 0.5 



+28 

-12 ± 

2 


-2.9 ± 0.5 

1993 February 

1993 March 

-20 

-10 ± 

2 

29 

-4.6 ± 0.9 



-10 

-6 ± 

4 


-2.9 ± 1.9 



+ 14 

-20 ± 

6 


-9.5 ± 2.8 

1993 March 

1993 April 

-20 

-10 ± 

2 

28 

-4.8 ± 1.0 



+ 14 

-24 ± 

4 


— 11.1 ± 2.0 



+24 

-24 ± 

4 


-11.1 ± 2.0 


630 


SCIENCE 


VOL. 265 • 29 JULY 1994 



influences and, among other things, to de- 
termine the origin of the hemispherical 
asymmetries in the temperature field. Other 
questions might be answered by a combina- 
tion of dynamical models and continued 
observations through at least another sea- 
sonal cycle, particularly in view of the po- 
tential for improved data quality. Several 
such questions are associated with zonal 
wave structure. Our results indicate upward 
propagation of tropospheric waves to the 
stratosphere in late 1988 (Figs. 7 and 8), but 
other zonal waves had no apparent influence 
on the stratosphere. Did the 1988 waves 
propagate upward solely because they coin- 
cided with a weak vertical shear over the 
NEB? Do the regular wave structures ob- 
served by Deming et al. (30) extend to the 
temperature field alone? There is clearly 
some indication of slowly thermal waves in 
the NEB and SEB; what determines the 
dominant wavelengths and phase speeds of 
these waves? Understanding wave behavior 
at lower latitudes may be very significant in 
studying the Galileo probe results for the 
temperature profile near 7°N latitude. 

The vast improvements in efficiency of 
recent technological developments in mid- 
dle-infrared array detectors compared to ras- 
ter-scanning and middle-infrared cameras, 
should provide greater coverage of Jupiter 
over a broader spectral range, characterized 
by higher signal-to-noise ratios with much 
better prospects for absolute flux calibration. 
Although angular resolution in the 18- pm 
region is diffraction-limited at the IRTF, 
observations with larger telescopes or appli- 
cation of image restoration techniques may 
improve the effective spatial resolution, en- 
abling us to study smaller features in the size 
regime where numerous visible and near 
infrared features illustrate the full variety of 
Jupiter’s diverse mesoscale weather. 
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